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Abstract 

This paper is concerned with the inverse problem of scattering of time-harmonic acoustic waves 
from a penetrable and buried obstacles. By introducing a related transmission scattering problem, a 
Newton iteration method is proposed to simultaneously reconstruct both the penetrable interface and 
the buried obstacle inside from far-field data. A main feature of our method is that we do not need to 
know the type of boundary conditions on the buried obstacle. In particular, the boundary condition 
on the buried obstacle can also be determined simultaneously by the method. Finally, numerical 
examples using multi-frequency data are carried out to illustrate the effectiveness of our method. 

1 Introduction 

In this paper, we consider the problem of scattering of time-harmonic acoustic waves from a penetrable 
obstacle with a buried impenetrable obstacle inside. Such problems occur in many applications such as 
radar, remote sensing, geophysics and nondestructive testing. 

Let Sq e C 2 denote a simple closed smooth curve in K 2 . Then R 2 is divided into the unbounded part 
Q.q and the bounded part Q. by Sq. We assume that Qo is filled with a homogeneous medium with the 
constant refractive index 1 and Q is filled with another homogeneous medium with the constant refractive 
index n\ > 0. Assume further that Q.2 cc Q. is an impenetrable obstacle with a C 2 boundary S \ and let 

Qi = o.\a 2 - 

Suppose an incident plane wave u'(x,d) = e ik ° xd , where ko = o)/c > is the wave number with oj 
and c being the wave frequency and speed in Qq an d d e S 1 is the incident direction, is incident on the 
penetrable obstacle Q. from the unbounded domain Qo. Then the total field u = u' + u s , which is the 
sum of the incident wave u 1 and the scattered wave u s in Q.q, satisfies the following Helmholtz equations 
together with the transmission condition on the interface S o and the boundary condition on the boundary 
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Figure 1: Scattering by a penetrable and a buried obstacles 
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where ki > is the wave number in Q.\ given by k 2 { = kxni, Aq is a positive constant depending on the 
property of the medium in Q. and Q.o, v is the exterior normal vector on So and +/- denote the limit from 
the exterior and interior of the boundary, respectively. The boundary condition SS{u) = represents the 
Dirichlet or Neumann boundary condition depending on the physical property of the obstacle ^2- It is 
well-known that the Sommerfield radiation condition (11.51) implies the following asymptotic behavior of 
the scattered field u s (see 0): 



where x = x/r is the observation direction and Uoo is called the far field patten of the scattered field u s . 

By using a variational or integral equation method it has been proved in ll22l |23| that the scattering 
problem dl.lMl.5b has a unique solution. It was further proved in E21 that So, Si and the boundary 
condition on Si can be uniquely determined by the far-field pattern iico(x,d,ko) for all Jc, d e S and 
a fixed ko if Aq ^ 1 is known. In this paper, we are interested in the numerical reconstruction of the 
interface So and the buried impenetrable obstacle Si together with the boundary condition on Si from 
the far-field pattern Moo(x, d, ko), given Aq > and the refractive index ni > 0. 




r = \x\ — > oo, 
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Many reconstruction algorithms have been developed for the numerical reconstruction of the inter- 
face Sq in the case when the buried obstacle flo = ( see > e -g-> the iteration methods in (HOT) and 
the singular sources method in Coll . the linear sampling method in [25] and the factorization method in 
lTT4l[T5l for the case when Aq = 1). In the case when the interface 5o is known and Aq - 1, a reciprocity 
gap functional method was introduced in ||3]|4j|6j|7l to reconstruct the buried obstacle Q.q from near field 
Cauchy data, whilst in [29] an iteration method was proposed to recover the buried sound-soft obstacle 
from the far field data. Recently, for the case when Qo ± 0, Cakoni et al. [2] introduced a multistep reci- 
procity gap functional method for reconstructing the interface S q and the buried obstacle flo recursively 
from near field Cauchy data under the condition that the boundary condition on S i or Q.q is known in 
advance, whilst in [30] the factorization method was extended to reconstruct the interface Sq from the 
far field data, without knowing the buried obstacle Qo- 

In this paper we develop a Newton-type iterative method to simultaneously reconstruct the interface 
5 o and the impenetrable obstacle S i together with the boundary condition on S i from the far-field pattern 
Uco(-, d, ko) for a finite number of incident directions d = dj, j = 1, P and a finite number of frequencies 
or wavenumbers ko = koi, I = 1, Q, given Aq > and the refractive index n\ > 0. A main feature of 
our method is that the boundary condition on the buried obstacle S i does not need to know in advance; 
in fact, the boundary condition on S i can be determined by the iterative method. This is different from 
the previous iteration methods which require to know the boundary condition on the obstacle to be 
reconstructed as a priori information. To simultaneously reconstruct Sq,S i and the boundary condition 
on S i, we first introduce a related transmission scattering problem by replacing the impenetrable obstacle 
Qo with a penetrable one with an arbitrarily chosen wavenumber &2 and having a transmission condition 
on S\ with the unknown transmission constant A\. We then use a Newton method to simultaneously 
reconstruct S q, S i and A\ . The boundary condition on S i can thus be determined in terms of the value of 
the reconstructed A\. See Section|2]for detailed discussion. 

Over the last 20 years, many work has been done related to Newton-type methods. For example, |[T6l 
gives the first rigorous characterization of Frechet derivatives of the far-field operator corresponding to 
an impenetrable, sound-soft obstacle embedded in a homogeneous medium by employing a variational 
method, and 11271 obtained the same result via the integral equation method. For the other cases such as 
impedance and transmission problems, we refer to ifTOlfTTI . It should be pointed out that the convergence 
of Newton-type methods for inverse scattering is far from complete. For partial results, the reader is 
referred to |[TT1 [T2l [T3l l28l . In addition, much effort has been made to reduce the computation cost of 
the Newton-type methods, such as the hybrid method lPT9l l20l and the methods of nonlinear integral 
equations lfTl[T8l. 

This paper is organized as follows. In Section [2j a related transmission scattering problem is pre- 
sented. In Section [3]the integral equation method is proposed for solving the related transmission prob- 
lem, and some quadrature rules are discussed briefly. The Newton method is presented in Section [4] and 
numerical examples using multi-frequency data are given in Section[5]to illustrate the effectiveness of our 
method. In the appendix, we give characterizations of the Frechet derivatives of the far field operators 
for the related transmission scattering problem in Section |2j employing the variational method. 

2 A related transmission scattering problem 

As discussed in the previous section, to simultaneously reconstruct Sq,S\ and the boundary condition on 
S i , we consider the following transmission scattering problem which is related to our original scattering 
problem: find the total wave u := u' + u s (which is the sum of the incident field u' and the scattered field 
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u s in Qq) such that 



Am + IcqU = in Qq 
Am + /cjM = in Qi 
Am + k\u = in Q2 



(2.6) 
(2.7) 
(2.8) 



M + - M 



M + — M_, 




UU+ UU- 

-7— - ^1-7- on Si 



(2.10) 



(2.9) 



where m', &o> &i, ^0 are as defined in SectionQ] k 2 is an arbitrarily chosen positive constant, A\ is a constant 
to be determined, and u s satisfies the Sommerfeld radiation condition (I1.5I ). This problem is obtained 
from the scattering problem (|1.1| >- (|1.5I > by regarding the impenetrable obstacle Qo as a penetrable one 
having an artificial wavenumber k 2 and a transmission condition on its boundary S 1 with the unknown 
transmission constant A\. 

We observe that the Dirichlet and Neumann boundary conditions can be considered as the limiting 
cases A\ — > 00 and A\ — > 0, respectively, of the transmission condition (12.10I ). Based on this observation, 
we will develop a Newton method to simultaneously reconstruct So, S 1 and A\ and then use the value of 
the reconstructed Ai to determine the type of boundary conditions on the boundary S 1 of the impenetrable 
obstacle Q.q. Precisely, if the value of the reconstructed A\ is very close to zero, then the boundary 
condition on S 1 will be regarded as a Neumann condition. If the value of the reconstructed A\ is far away 
from zero, we define t\ := l/A\, U := A\u\q 2 and replace (12.81 ) and (12- 10b . respectively, by 



We then use the Newton method to simultaneously reconstruct So, Si and t\. If the value of the re- 
constructed T\ is very close to zero, then the boundary condition on S 1 will be regarded as a Dirichlet 
condition. 

Remark 2.1. During the preparation of this paper we were notified by the authors of lETj that they 
have recently also studied the asymptotic behavior of the solution to the scattering problem by a lossy 
inhomogeneous medium buried in a layered homogeneous medium as the medium mass density is very 
large in the context of cloaking or invisibility; precisely, they proved mathematically that the solution 
of the transmission problem (I2.6I) - (I2.10I ) and (11.51 ) converges to the solution of the scattering problem 
dl.lb - dl.5b with a Neumann boundary condition 33(u) := du/dv = on S 1 as A\ — > 0+ in the case when 
Im(^) > 0. 

Based on the above observation, we assume that ko, k\,k 2 , Aq are known and introduce the far field 
operator F which maps a set of admissible boundaries So and Si and the constant A\ into the far field 
pattern of the scattered solution to the corresponding transmission scattering problem (I2.6b - (I2. 10b and 
(1 1.5b . We also introduce the far field operator F which maps a set of admissible boundaries So and S 1 and 
the constant t\ into the far field pattern of the solution to the transmission scattering problem (I2.6b - (l2.7b . 
d2!9l) , (|2TTTT )- (|2~T21 ) and dT31) . We will solve the equation 



M + = T\ U-, —— - 



AU + kjU = in Q 2 

du+ dU- 



(2.11) 



(2.12) 



F(Sq,S\,A\) = u, 



00 



(2.13) 
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for Sq,Si, M; as discussed above, for the case when Ai is very large, we solve the equation 

F(S ,S 1 ,t 1 ) = u oo (2.14) 

for So,S i, Ti. Here, the value of A\ or ri reveal the type of boundary conditions onSi and Moo is the far 
field pattern of the original scattering problem dl.lMl.5l l. 

Following the idea of |[T6l [lOl . we can prove that if So and Si are C 2 boundary then F (or F) is 
Frechet differentiable with respect to So, S i and A\ (or ri) for the case when A\ € (0, oo) (or t € (0, oo)) 
(see Appendix for details). Motivated by this we propose a Newton iteration method to recover the 
boundaries So, Si and the boundary condition on S i in Section @] 



3 The integral equation method for the direct problem 

For the Newton method, we need to compute the Frechet derivative of F or F which is characterized 
by the solution of the problem d2.6M2.10l > and d 1 -5b with certain inhomogeneous boundary data (see 
Appendix). For convenience we consider the following more general transmission problem: find u s € 
C 2 (Q ) fl C 1,o "(n ) and u € C 2 (Qi) f| C l ' a (a{) f| C 2 (H 2 ) D C^fih) with < cr < 1 such that 

(3.1) 
(3.2) 
(3.3) 

du s , du- 

u + -u-=fu — A Q — = f 2 on S (3.4) 

dv dv 

du + dii- 

u + -u-- f 3 , — M — = / 4 on S i (3.5) 

dv dv 

and u s 'satisfies the radiation condition (fT31). where (fufij^fa) 6 X := C 1 ' cr (S )xC°' cr (So)xC 1 ' cr (S i)x 
C^Si). 

Remark 3.1. For our scattering problem d23 M2l0T) and (TOl . (/i , / 2 , / 3 , / 4 ) = (-m ! ', -du'/dv, 0, 0). 

The transmission problem d3.1| )- d331 ) and (11.51 ) has been studied in B41 . employing the integral 
equation method. In this paper we will also use the integral equation method to deal with the inverse 
problem. To this end, for / = 0, 1, 2 let Q>i{x, y) denote the fundamental solution of the Helmholtz equation 
A + kj in K 2 and, for i, j = 0, 1 define the boundary integral operators Syj, Kyi, KT. p Tyi by 
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(%^)(jc)= f Oi(x,y)^(y)ds(y) xeS,- 

l< 

_ 5v(y) 



(KwMx) = f m ' { ^\ {y)ds{y) xeSi 
JSj ' " " 



T C dO/(x,y) 

{Kj jl xlr){x)= J) \ yJ ifr(y)ds(y) xeSj 

1 JSj M*) 

(Tijif)(x) = f m ^\ (y)ds(y) xeSi 
dv{x) J s , dv(y) 
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We also define the far field operators Soo and Koa corresponding to Sjjo, Kjjo by 



Jx/4 



o X 



- ik ^-y^(y)ds{y) 



(K^fa(x) = 



e ix/4 r 

^Snk Jso 



Q e -ik x-y 

dv(y) 



&(y)ds(y) 



X € 



X € 



where S 1 is the unit sphere in R 2 . The reader is referred to [5 ] for mapping properties of these operators 
in the spaces of continuous and Holder continuous functions. 

Following 11241 . we seek a solution to the transmission problem (I3.1I )- (I3.5I ) and (11.51) in the following 
form: 



u s (x) 
u(x) 



u{x) 



X, 



&bo(x ,y) } + ^ Q(xy)lff2(y) 
dv(y) 

d®i(x,y) 



ds(y) 



X € Q, 







dv(y) 

3 d^(x,y) 
d® 2 (x,y) 



fa + ®i(x,y)fa(y) 



ds(x) 



foiy) + ®i(x,y)fa(y) 



dv(y) 



fa + ®2(x,y)fa(y) 



ds(x) 



ds(y) x e Q.\ 
x e D.2 



with (fa, fa, fa, fa) £ X. From the jump relations of the potential operators, it is known that the poten- 
tials u s and u defined above solve the problem (I3.1I) - (I3.5I ) and (11.51) if *P = (fa, fa, fa, fa) T solves the 
operator equation 



(/ + A)¥ = R 



(3.6) 



where R = (jiof\, -^0/2,^1/3, -^\f\) T , Hi - 2/(1 + Aj), i=0, 1, / is an identity operator and A : X — > X 
is given by 



^o(^o^ooo - ^ooi) 

-^0^0(^000 - ^ooi) -^o(^nnn ~~ ^o-^nni ) 



^o(^ooo - S 001) 
l ooi^ 



"000 
HiS 101 



101 



-j"o^i^oii -/A^on 

^o^o^i^on fioAoK on 

j"i(^i^in - ^112) l*i(S 111 - S 112) 

-[i\Ai(T\i\ - Ti\%) -n\(K\ n - A\K^ n ) , 



It was proved in |[24l that the system ( 13.61 ) of integral equations is uniquely solvable in X (see Theorem 
2.4 in [24] for details). Further, by the asymptotic behavior of the Green function 0, we have 



WooC*) = /loCSoo^lX*) + (KoofaXx), X. € S , 



(3.7) 



where (fa, fa) € C l,cr (So) X C°'°~(So) are the first two components of the unique solution to ( I3.6I ). 

To generate the far field data, we use the Nystrom method to solve the system (13.61 ) of integral 
equations. We assume that the boundaries S j,j = 0, 1, can be parameterized as x(t) = (x\(t),X2(t)) for 
< t < 2n, where x(t) is an analytic and 27r-periodic function. Then, by the change of variables the 
integral operators S iji, K^, Kf.j, i,j = 0,1,1 = 0,1,2, and Tjji (i,j = 0,1,1 = 0, 1,2, i + j) can be 
reduced to the form 



ijl 

-•271 



J |&i (t, t) In {Asin 2 -y- ) + h(t, r) 



(f)(T)dT 
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and Tjji (i = j) can be reduced to the form 



1 C r-t C \ I t — t\ 

— cot — 0'( T )dr + J o [*i(f, t) In |4im 2 — J + fc 2 (f, t) 

where t), ki(t, t) are analytic and 27r-periodic functions. 

In this paper we use the following quadrature rule: for < t < 2n, 



0(r)dr 



2n-l 



/^4sin 2 — }f(r)dT * Y^Rf{t)fitj) 
1 r 2?r x - f 2,1-1 

where, for / = 0, . . , , 2n - 1, 

# (n) (f) : = _^yicos[m(f-^)]-^cos[n(f-/ n) )] 
^ n ^—i m n l 3 

m=l 

1 " _1 1 

r (n) (0 := — Vffl cos[m(f - t (n) )] - - cos[n{t - t {n) )] 
J n z—L J 1 J 



m=l 



After finding the solution *F = (\fj\, if/2, fa, i/M) r to the system of integral equations, (13.61 ). we can use 
(IP) to compute the far field pattern (cf. 151 [1711291). 



4 The inverse problem 

In this section, we propose a Newton iteration method for the inverse problem. We assume that S o and 

5 i are starlike surfaces which can be parameterized by jq{6) and j\{9), respectively: 

7O (0) = (O,O) r + r o (0)(cos(0),sin(6O) r 
yi(0) = (a u a 2 ) T + r l (e)(cos(G),sm(6)) T 

where 9 € [0, 2tt] and (0,0) and (01,02) are their centers, respectively. We will consider the noised 
perturbation of the far field data « M in the sense that ||m£, - Mooll^^i) < <5||woollz2(si)> where 8 > is 
called the noise ratio. Thus, given w£,(-, di, ka), i - 1, ...,P for a fixed &o, we rewrite ( 12.131 ) and ( 12.141 ), 
respectively, as 

Fi<y ,7i,M) * < f , /=1,-..,P (4.1) 
^i(7o,n.ri) * <,-, i = l,...,P (4-2) 

where . = w£,0>^i>^b) an(1 ^i*-^; represent the far field operators (14.11) . (14.21 ) corresponding to the 
incident plane wave with direction dj. 
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We now derive the Newton method for (14.11) . The method is similar for (14.2b . Assume that, after the 
ra-th iteration, we have the approximate data y^\ y\ ' \ A^ with A^ ± 0, where 



yf\8) = (0,0)' +r^(cos(0),sin(0))' 
yf\Q) = (elf , aff + r^\9)(cos(9), sin(0)) 

Then the linearized form of (14.11) becomes 
l 



v El 

Zj fa,, 



dAi 



(y^^'J?) 



(4.3) 



for i = 1, ...,P, where Ayj m) (0) - yj'" +1) (<9) - yj m) (0), / = 0, 1, and A^ m) = ^ m+1) - A { ™ ] are the updates 



to be determined. The linearized equation (14.31 ) is ill-posed, so a regularization strategy is needed. Fol- 
lowing 11131 . we use the Levenberg-Marquardt algorithm, that is, we seek a solution (Ay|, m) , Ay{ m) , AA { ™ } ) 
with AA^ e R and 



Ay ( ™\G) = (0, Of + A^ n) {6){cos{8), sin{9)f 
Ayf\G) = (Aaf\Aaff + Ar^ (9)(cos(6), sm(6)) T , 

where Ar ( ™\ Ar{ w) e # 4 (0, 2tt), * > 0, and Aa[ m \ Aaf } e R, such that (AyJ* 3 , Ayj m) , Ai ( " ,} ) is a solution 
to the minimization problem: 



Zj Zj 5 r/ 



;=i /=i 



(Ay (m) ) + ^ 



3F; 



(A/tj ) + / , ,(y () ,yj ,Aj ) - 



^\\Arf% s(Sl) + ^\Aa^\ 2 + {AA^ 



i=i 



i=i 



2 

(4.4) 



Here, /? € K. + is chosen such that 
dFi 



Zj Zj ^v, 



V i=l /= 



(Ay«) + 



3F; 



dAy 

Zllz7/ ( m ) ( m ) l( m )\ ' 
||F,-(y^ ',y\ ',A\ >)-u 



(A/tj ) + J r,(y ,yj ) - u^j 



( p 



\i=l 



2 

lHs 1 ) 



1 

\2 



LHS 1 )) 

(4.5) 



with a given constant p < 1. After obtaining (Ayi , Ay^ , ), the new approximation to (yo,yi,^i) 
will be (yg" +1) , y< m+1) , ^' n+1) ) = (y^ + AyJ' 1 ', yf ] + Ayf\ A ( ™ } + AA { ™ } ). In the numerical computation, 
we can use the bisection algorithm to determine (see |[T3l ). Note that is unique if certain assumptions 
are satisfied (see [8 ]). In this paper we choose f5 by a posterior parameter choice rule. For other strategies 
of solving (14.11 ) we refer to lfT3l . 

Further, we define the relative error by 



II 17 / (m) (m) ilmk 

i p \\ F >(y >7i ) 



Oh) 
ko 



lHs 1 ) 



(4.6) 



LHS 1 ) 
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Then the iteration is stopped if Err, < t5, where t > 1 is a given constant. Here, we use the subscript 
to emphasize its dependance on the wave number /co- 
in the numerical experiments, for fixed d, ko we use ^(i,, d, ko), i = 1, to denote the mea- 
sured far field data, where f,- = 2n(i - l)/n. We then seek the approximate solution y ,7j , A™ with 
r m) , r ( j m) € R M c H\0, In) (s > 0) and A^ e R, where 

M 

R M := \r € //' s (0, In) \ r(9) = a + ^ [a/ cos(/<9) + a/ +M sin(/<9)], a/ € R, I = 0, 2M}. 

i=i 

Similarly as in [T3fl, we use the following: 

H/fesi) * — £l/(*«)l 2 , /eC(S') 

M 

\\g\\ 2 HS(Sl) := 27ra + ^2][(l +/ 2 )>? + ^ +M )] 

/=i 

where g(ff) = ao + [f / cos(W) + ai + M sin(/0)]. Then (I4.4I ). (14.5b . (14.61 ) can be approximated, and the 
updated functions are computed by determining their coefficients. For the determination of the coeffi- 
cients the reader is referred to [13], and for the computation of the Frechet derivative of the far field 
operator we refer to lfTTl[T2l[T6l . 

Remark 4.1. In Section |2j we observed that if A\ — > then the solution u of the transmission problem 
(I2.6I )- (I2.10I ) and (11.51 ) converges to the solution of the scattering problem (ll.ll )- (ll.5l ) with a Neumann 
boundary condition 38{u) := du/dv = on Si and that if Ai — > oo (or t\ —> 0) then the solution u of 
the transmission problem (12.6I) - (12.10I) and (11.51 ) converges to the solution of the scattering problem (11.11 )- 
(11.51) with a Dirichlet boundary condition 38(u) := u = on Si. Thus, the following fact is expected: 1) 
when Woo is the far field pattern corresponding to a sound-hard obstacle Q.2 (i-e., a Neumann boundary 
condition on S i) then A^ given by the Newton iteration method for (14.11 ) will be very close to zero, and 
2) when Woo is the far field pattern corresponding to a sound-soft obstacle Q.2 a Dirichlet boundary 
condition on S \ ) then t ™ given by the Newton iteration method for (14.21) will be very close to zero. 
However, there is no rigorous proof of this convergence result even for the case of an obstacle in the 
homogeneous background medium. 

Remark 4.2. Since we do not know the boundary condition on S i , we choose a constant A. Let A j m) be 
the approximation of A\ after the m-th iteration and T^ m) = \/A ( ™\ Then, if |A^ m) | <= A, we solve (|4~TT) to 
get (r m+1) ,rf +1 U5 m+1) ) and let t ( ™ +1) = l/A ( ™ +1) . On the other hand, if |^ m) | > A then we solve g2) 
to get {Jq , 7j , Tj m+1) ) and let /lj" !+1 ^ = 1/tj I It should be noted that, in numerical computation, 
A^ may be negative. This does not affect the iteration process, and A^ can still be used to determine 
the type of boundary conditions on S i . 

Remark 4.3. In our iteration algorithm, multiple frequency data are considered in order to get better 
reconstructions. General speaking, the low frequency data are used to get the main profiles of the bound- 
aries and the higher frequency data are used to get a refined reconstruction. 

Our Newton iteration algorithm is given in Algorithm |4J] below. 

Algorithm 4.1. Given Aq, di, koj), i - 1, P, j = 1, Q, where &oi < ^02 < • • • < k^p. 
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1 ) Choose initial guesses 7 \ y\ > ^\ far jq, j\ , A\, where A^ + 0, and let Ir be the maximal iteration 

number far each fixed frequency. Set j = 0, = 1 and go to Step 2 ). 

2) Set j - j + 1. If j > Q, then stop the iteration; otherwise, set Icq - Icqj, m = and go to Step 3). 

3) If Err^ < tS or m - Ir, then set (y^,y^,A^ ) = (Jq, 7^"\ A j m \ Tj" ) and return to Step 2; 

otherwise, go to Step 4). 

4) Using the strategy in Remark \4~2\ we get the new approximate data (yn m+1 \ Ji \ A ™ + , t^™ + ). 

5) Set m = m + 1, and return to Step 3). 

Remark 4.4. Let A^ nd) be the final approximation of A\. If l/l^"^! >> 1, then it is concluded that S i is 
a sound-soft obstacle, and if |/l ( .™ rf) | « 1, we then conclude that S \ is a sound-hard obstacle. 



5 Numerical examples 

In this section, we present several numerical experiments to demonstrate the effectiveness of our algo- 
rithm. In all numerical experiments we make the following assumptions. 

1) We choose n x - 0.64,i = 1.2, s - 1.6 (which is suggested in US), M=25, p = 0.8, t = 1.5 

and A = 1 for iteration. Since ki is a parameter at our disposal, we just choose kj = k\ for each 
frequency. 

2) To generate the synthetic data, we use the integral equation method in Section [3] to solve the direct 

problem. In order to avoid inverse crimes, we choose n = 128 for generating the synthetic data 
and n = 64 for solving the direct problem at each iteration. 

3) We use the full far-field data with 64 measurement points. The noisy data u^ is obtained as w^, = 

Moo + (^llMooll^gij/II^H^gi), where £ is a random number with Re(£), Im(£) e A^(0, 1). 

4) The initial shapes of So, S i are chosen to be circles. 

5) For the incident plane waves, we choose d[ = (cos(6i), sin{9i)), i = 1, ...,P, where P is the number of 

incident plane waves and = 2n(i - 1)1 P. 

6) In each figure, we use solid line dotted line and dashed line '- -' to represent the actual curve, 

the initial guess of the curve and the reconstructed curve, respectively. 

7) The shapes of the actual boundaries are given in Table [Q 



Type Parametrization 



Circle ro(cos t, sin t), t e [0, 2n] 

Apple shaped [(0.5 + 0.4 cos t + 0.1 sin(20)/(l + 0.7 cos f)](cos t, sin t), t e [0, 27r] 

Kite shaped (cos t + 0.65 cos(2f) - 0.65, 1.5 sin 0, t e [0, 2n] 

Rounded square (3/2)(cos 3 t + cos t, sin 3 t + sin t), t e [0, 27r] 

Rounded triangle (2 + 0.3 cos(3f))(cos t, sin t), t € [0, 2n] 



Table 1: Parametrization of the boundaries So and Si. 
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Example 1. In this example, we consider a sound-soft, apple-shaped obstacle embedded in a rounded 
triangle-shaped penetrable obstacle. Suppose the far-field data are given for this case which corresponds 
to a Dirichlet boundary condition on the boundary S i of the embedded obstacle. However, the boundary 
condition on S i is assumed unknown, and we will use our method to reconstruct both obstacles and 
determine the type of boundary conditions on S i . For the initial guesses y®\ y^ of 5 o, S i , respectively, 
we choose the radii rjj^ - 2.4, rf^ = 0.5, respectively, and the two centers to be the origin (0,0). 
By trials, choose /l (0) = 10. Here, the center of y^ is updated at each iteration. In Figure |2j the 
reconstruction results are given by using the exact data with a fixed wave number Icq = 2. Figure |2(a) 



|2(b)[ |2(c)| and |2(d)| show the reconstructed results with one, two, three and four incident plane waves, 
respectively. It is observed that using more incident plane waves yields a better reconstruction result. 
Figures |3] and |4] give the reconstruction results by using the 3% noisy far-field data with one and four 
incident plane waves, respectively. Here, we use multiple frequency data with Icq = 2, 4, 6, 8. From 
Figures [3] and @] it is seen that multiple frequency data yield excellent reconstruction results. Finally, 
since \A\ \ » 1 for each case, we can conclude that 5 1 is a sound-soft boundary. 

Example 2. In this example, we consider the problem in Example 1 again. Here, the far-field data are 
given again for the case of an embedded, sound-soft obstacle, and S i is also assumed to be a sound-soft 
boundary. We will use the Newton iteration method for the inverse problem corresponding to the original 
scattering problem dl.lMl.5l l with 3§{u) := u to reconstruct the two obstacles. The initial guesses of 
yo,y\,M) are the same as in Example 1, and the center of Y\ is updated at each iteration. In Figure [51 
we present the reconstruction results by using the 3% noisy far-field data with one incident direction and 
multiple frequencies &o = 2, 4, 6, 8. From Figures |3j and |5jit can be seen that the reconstruction obtained 
by both methods is almost the same. 

Example 3. In this example, we consider a sound-hard, kite-shaped obstacle embedded in a rounded 
square-shaped penetrable obstacle. We assume that the far-field data are given for this case which corre- 
sponds to a Neumann boundary condition on the boundary S i of the embedded obstacle. However, the 
boundary condition on S i is assumed unknown, and we will use our method to reconstruct both obstacles 
and determine the type of boundary conditions on 5 1. For initial guesses of S o, S i, we choose the radii 
ri = 2.8, r \ — 1 and the two centers to be (0, 0) and (0.5, 0), respectively. By trial, we choose = 1. 
Since So and S i are very close, it is more difficult to reconstruct both obstacles. In order to get a stable 
reconstruction, we use at least three incident plane waves. Moreover, the center of is updated only 
at the first five iterations. Figure [6] shows the reconstruction results using the exact far-field data with a 
fixed wave number Icq - 1. Figures [6(a)[ |6(b)| |6(c)| and |6(d)| give the reconstruction results with three, 
four, eight and sixteen incident plane waves, respectively. It is observed that using more incident plane 
waves does not yield a better reconstruction even for the exact data. In Figure [7] we present the recon- 
struction results using the 3% noisy far-field data with three incident directions and multiple frequencies 
ko - 1, 3, 5, 7. Compared with Figure[6l using multiple frequency data gives a much better reconstruction 
even for the noisy data. Since \A { * r) \ « 1 for each case, we conclude that S i is a sound-hard boundary. 

Example 4. In this example, we consider the same problem as in Example 3. Here, the far-field 
data are given again for the case of an embedded, sound-hard obstacle, and S i is also assumed to be a 
sound-hard boundary. We will use the Newton iteration method for the inverse problem corresponding 
to the original scattering problem dl.lMl.5l > with £%{u) :— du/dv to reconstruct the two obstacles. The 
initial guesses of yo, yi,A$ are the same as in Example 3, and the center of y^ m) is updated only at the 
first five iterations. In Figured! we present the reconstruction results by using the 3% noisy far-field data 
with three incident direction and multiple frequencies ko = 2,4,6,8. Figures |7] and [8] indicate that the 
reconstruction obtained by both methods is almost the same. 
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(a) no noise, A ( ' r) = -502.6, Err (lr> = 2.151 x 1(T 3 , Ir = 25 (b) no noise, A ( ' r) = 1004, Err (lr) = 2.246 x 10~ 3 , Ir = 25 




-3-2-10123 -3-2-10123 



(c) no noise, A { ' r) = -9971, Err <Ir) = 1.861 x 10~ 3 , Ir = 25 (d) no noise, A ( ' r) = 2796, Err ilr) = 1.975 x 10~\ Ir = 25 

Figure 2: An apple-shaped, sound-soft obstacle embedded in a rounded triangle-shaped penetrable ob- 
stacle. The boundary condition on S i is assumed unknown, and the exact data with a fixed wave number 
ko = 2 are used, (a), (b), (c) and (d) give the reconstructed results with using one, two, three and four 
incident plane waves, respectively. 
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(a) 3% noise, k = 2, A ( ' r) = -275 .4, Err (Ir) = 0.04053, (b) 3% noise, k = 4, /if = -515, Err {lr) = 0.03803, /;• = 6 
Ir = 12 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 



(c) 3% noise, &o = 6, A ( ' r) = 736.2, Err Ur) = 0.04316, /r = 2 (d) 3% noise, *o = 8, % r) = -895.7, Err ( " } = 0.03651, 

/r = 4 

Figure 3: An apple-shaped, sound-soft obstacle embedded in a rounded triangle-shaped penetrable ob- 
stacle. The boundary condition on S i is assumed unknown, and noisy multi-frequency data with one 
incident plane wave are used, (a), (b), (c) and (d) show the reconstructions at &o = 2, 4, 6, 8, respectively. 
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-2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 



(a) 3% noise, k = 2, Af' } = -2591, Err Ur) = 0.03686, (b) 3% noise, k = 4, A { ' r) = 963.5, Err {lr) = 0.03987, Ir = 4 
Ir =12 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 



(c) 3% noise, k = 6, A C / r) = -417.3, Err'" 1 = 0.04476, (d) 3% noise, fc = 8, A { ' r) = 5715, Err {lr) = 0.04139, /r = 2 
/r= 2 

Figure 4: An apple-shaped, sound-soft obstacle embedded in a rounded triangle-shaped penetrable ob- 
stacle. The boundary condition on 5 1 is assumed unknown, and noisy multi-frequency data with four 
incident plane waves are used, (a), (b), (c) and (d) show the reconstructions at Icq = 2, 4, 6, 8, respectively. 
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(a) 3% noise, Err (lr) = 0.04301, Ir = 12 



(b) 3% noise, Err (,r> = 0.04308, Ir = 5 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 



(c) 3% noise, Err (lr) = 0.03717, Ir = 3 (d) 3% noise, Er/' r) = 0.03687, Ir = 2 

Figure 5: An apple-shaped, sound-soft obstacle embedded in a rounded triangle-shaped penetrable obsta- 
cle. The boundary condition on S \ is assumed known, and noisy multi-frequency data with one incident 
plane wave are used, (a), (b), (c) and (d) show the reconstructions at ko = 2,4, 6, 8, respectively. 
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(a) no noise, A { ' r) = 4.186x 1(T 4 , Err {lr) = 0.001296, Ir = 30 (b) no noise, A ( ' r) = 2.097 x 10~ 5 , Err {,r) = 0.01086, Ir = 30 




(c) no noise, A { ' r) = 6.035x 10~ 4 , Err ( ' r > = 0.001 196, Ir = 30 (d) no noise, Af ] = 4.963x 10~ 4 , Err {,r) = 0.001212, Ir = 30 

Figure 6: An kite-shaped, sound-hard obstacle embedded in a rounded square-shaped penetrable obsta- 
cle. The boundary condition on S i is assumed unknown, and the exact data with a fixed wave number 
ko = 1 are used, (a), (b), (c) and (d) show the reconstruction results using three, four, eight and sixteen 
incident plane waves, respectively. 
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(a) 3% noise, k = 1, A { ' r) = -1.085 x 1(T 3 , Err ( "> = (b) 3% noise, k = 3, /tf' = -3.155 x 1(T 3 , Err i,r) = 
0.04189, Ir = 14 0.04154, Ir = 9 




-3-2-10123 -3-2-10123 



(c) 3% noise, k = 5, ^ ( , /r) = 4.916x Err i,r) = 0.04464, (d) 3% noise, k = 7, A ( / r) = -6.377 x 10" 3 , Err*" 1 = 
Ir = 6 0.04146, Zr = 2 

Figure 7: An kite-shaped, sound-hard obstacle embedded in a rounded square-shaped penetrable ob- 
stacle. The boundary condition on 5 1 is assumed unknown, and noisy multi-frequency data with three 
incident plane waves are used, (a), (b), (c) and (d) show the reconstructions at Icq = 1, 3, 5, 7, respectively. 
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-3-2-10123 -3-2-10123 



(a) 3% noise, k = 1, Err (,r) = 0.04124, Ir = 13 (b) 3% noise, k Q = 3, Err {,r) = 0.03814, Ir = 9 




-3-2-10123 -3-2-10123 



(c) 3% noise, k = 5, = 0.03943, Ir = 6 (d) 3% noise, k Q = 7, Err*"'' = 0.04089, Ir = 2 

Figure 8: An kite-shaped, sound-hard obstacle embedded in a rounded square-shaped penetrable obsta- 
cle. The boundary condition on S i is assumed known, and noisy multi-frequency data with three incident 
plane waves are used, (a), (b), (c) and (d) show the reconstructions at Icq = 1, 3, 5, 7, respectively. 
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Appendix: Frechet derivatives 



In this appendix we study the Frechet derivative of the far field operators F and F when A\ > and 
T\ > 0, respectively. 

Let (ho,h,AA\) e C\Sq) x C\S i) x R be small perturbations of (So,Si,A\) (that is, 
||^|| c i (Si) , A/li « 1, i = 0, 1). Define (S ,•)/*, := {y e R 2 \y = x + hfa), x e S,}, i = 0, 1. The far 
field operator F is called Frechet differentiable at (So,S\,Ai) if there exists a linear bounded operator 
F'\(s , Sl M) ■ C\S ) x C l (S i) x R -» L 2 (S ') such that 

l|F((So)/!o,(S + A ^i) - ^(5 ,5i,/li) - F'| (SojSlj/il) (/i ,/ii, A^i)|| L 2 (S i } 

-o(||/2o|lc'(5o) + ll /z illci(s 1 ) + l A ^il) 

^"l(S .5i.'li) * s cai l e d the Frechet derivative of F at {Sq,S y,A\). Following the idea in ifTOl 1X61 . we have 
the following theorem to characterize the Frechet derivation of F when A\ > 0. 

Theorem A.l. Let So, S \ be C 2 and let u € H} oc (R 2 ) be the solution of the problem ri2.6D - rf270l ) and 47.51 ) 
with A\ > 0. 77je7i F is Frechet differentiable at(So,S\,A{) and the Frechet derivative F'\^ Q ,s iJO w g/verc 
fry ^'l(5o,Si^i)(^o» h\,AA\) = u'cc, where u'^ is the far-field pattern ofu' e H} (Qq) f] H l (Q.\) f] H (Cl?) 
which solves the problem dJ.il )-f T01 ) and rti.5D w/?/j 



/l = 


(du+ du-\ 


So' 








h -- 


- - A k 2 )(h Q )yu\ SQ 


+ Div So 


[(h ) v ((Vu + ) t 


- A (Vu. 


-)*)], 


h -- 


(du+ duA 
'- - {h) \-dv- ~ -dv-) 


Si' 








f4 -- 


- (k 2 - A\kl)(hi) v u\ Si 


+ Div Sl 


[(/ii)v((Vu + ), 


-Ai(Vu. 


-),)] + 



AAi du^ 



A\ dv 

where h v and h t denote the normal and tangential components of the vector field h. 

Remark A.2. From (9l it follows that u e H 2 oc (D.q) f] H 2 (Q.\) f] H 2 (Q.2), which guarantees the existence 
of u' . 

Proof. Our proof is based on the variational formulation. To this end, let Br be a fixed ball of radius 
R and centered at the original with R large enough such that Q.\ [j Cl 2 c Br and let L be the Dirichlet 
to Neumann (D-N) map defined on 3Br [5]. It is easy to see that the solution u satisfies that for all 
v e H\B R ), 



(aVu ■ Vv - k 2 uv)dx- < Lu 



JdB R \ dv 



va 



Lu j vdx 



(A.l) 



with 



a 



1 in Q.Q, 
A in Qi, 
AqAi in Q, 2 , 

-1/2 







Ank 2 



AqA\ k 2 

1/2/ 



in Q,q, 
in Qi, 
in Q.2, 



where < •, • > denotes the dual between FT ' (3Br) and H ' (3Br). For the property of L and the proof 
of the equivalence between the variational equation and the origin problem, we refer to ll23Tl where a 
similar problem is considered. 
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Denote by S(u,v) the left-hand side of (IA.1I ). Then, by the Fredholm alternative and the uniqueness 
of the origin problem, it can be proved that there exists a linear bounded operator B : H\B R ) — > H l (B R ) 
with a bounded inverse such that 



SQt,v) = (B(u),v), 



(A.2) 



where (•, •) denotes the inner product in H l (B R ). Thus, for any / € H l (B R ) there exists a unique function 
u f e H l (B R ) such that S(u f ,v) = (f,v) for all v € H\B R ). 

Let Kq,K\ be two domains with So c Kq cc B r \D. 2 , Si c K\ cc B r \Q.q and Kq n K\ = 0. 
For (ho,h\) e C 1 (So) x C (Si) there exists an extension of QiQ,h\) which is denoted again by QiQ,h\) 
such that ho,hi e C l (R 2 ), supp(/jo) c Kq, supp(/ji) c and Pilled) - CPi'llcks,-)' ' ~ ^' ^ ^ et 
/j(x) - ho(x) + h\(x) for x € K 2 and let y>(x) = je + /j(x). If p,- licks-)! i - 0, 1, are small enough, then 
ip : B R — > Br is a diffeomorphism (see ifTBl ") and ^ maps So to (So)/ 10 , S i to (Si)/,, , respectively. Denote 
by \p the inverse of <p and by J v , the Jacobians of ip, iff, respectively. 

Let U(h,AAi) t> e tne solution of the problem (I2.6I >- (I2.10I ) and (I1.5I) with (So,Si,Ai) replaced by 
((So)ho> (S > ^i + A/li). Then U(hMi) satisfies that for all v € H l (B R ), 



Ibr 



— 2 



r (du!_ . 

JdB R \ dv 



vdx 



(A.3) 



with 



1 in (Q )a, 
^o in (Oi)^, 
/l (/li + Aii) in (fl^)*. 



&o in (Q )/„ 
/l ^i in (Qi)ft, 
i (^i + A/li)^ in ( Q 2)ft, 



where (Q,-)fc = + h(x)\x € O,}, / = 0, 1,2. If A/li is small enough, then the unique solution of 
the problem (IA.3I) exists and can be extended outside of B R such that m^aai) £ ^hc^ 2 ^ w ^ tn "(M^i) 
satisfying ( 12.61) and w^Aii) _ " ! satisfying dl.51 ). 

Write U(h,Mi) ~ u (h,AA,) ° f- Then, by the change of variables, we obtain that for all v € H l (B R ), 



Ib r 



(aVSy ! Aii) 7^7^Vv - ^u^A^detiJ^dx 



+A AAi (yuJ hMi) J,j,J^Vv - klu( hMl) v)det(J v )dx- < Lu (hMl) ,v >= J - Lu'J vdx 

Denote by S^A/yO^A/h), v ) the left-hand side of the above equation. Then, for all v € H l (B R ) we have 

S {u - U(hM\)' V ) = $ (h,AAl)(U(hMlh v ) ~ S(U(h,AAl)' V ) 

= I [aVii{ hMi) {]$J\ det(7 y ) - /) Vv - k 2 (det(7 y ,) - l) S ( / 1jA/ll) v] dx 

J Br 

+AqAA\ J (VS^ A/ll) 7^7^Vv - ^"(/i.Aiov) det(/ ¥ ,)<ix 
Let w € H l (B R ) be the solution of the problem: find w e H l (B R ) such that 

5(w,v) = J~ [arVw r (//, + /[ - V • (W)) Vv + k 2 (V ■ h)uv] dx 

-AqAAi J (Vu T Vv-k 2 2 uv)dx, Vv € H l (B R ) 



(A.4) 
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(by the property of S(-, •), it is easy to prove that such a w exists). Then, and since \J#f£ det(7 v ) - I + 

(J h + J\ - (V • h)I)\ = o\\\h\^ -J, det(7 v ) - 1 = V • h + o(\\h\\ 2 cl -J and tends to u in 

(which can be proved similarly as in [ 10]), it is easy to prove that 



S(u {hMl) -u-w,v) / v 
sup — = {\\ h WcHB- R ) + M ) ■ (A.5) 

Then, by CO> we have [l/(\\h\\ c i m + |AA 1 |)](S ( / ki yu 1 ) - u - w) -> in H Y (B R ) as ||A|| c i (fljd + IMI tends 
to zero. This, together with the trace theorem, implies that [I/CPHcVbjj) + \^M\)](u(h,AAi) - « - w) — > 
in H lj2 (dB R ) as p|| c i( Bi ) + \AA\\ tends to zero. 

Now let u' = w - h ■ Vm. Then u' e H l (B R n Q ) f^H 1 ^) f}H l (Q. 2 ) since u e 
fl^(fio) # 2 (^i) ^ 2 (^2)- Note that 

Vu T [J h + J\ - (V • h)I]Vv = V • [(ft • Vv)Vm - (Vm • Vv)/j] + V{h ■ Vm) • Vv - (h • Vv)Am. 

Then for v € H 2 (B R ) the right-hand side of (1A.4I) becomes 

I [qtV(/z • Vm) • Vv + aV • ((h • Vv)Vm - (Vm • Wv)h) - a(h ■ Vv)Am + & 2 (V • /i)hv] dx 
Jb r 

-AqAAi J (Vm t Vv - k\uv)dx 

= I [arV(/z • Vm) • Vv - & 2 (/i • Vm)v + aV • ((ft • Vv)Vm - (Vm • Vv)ft) + & 2 V • (ftwv)l die 

-AqAAi I V • (Wuv)dx. 
Jn 2 

On the other hand, the left-hand side of (IA.4b is 5(w,v) = f„ (uVw • Vv - k 2 wv)dx- < Lw,v >. Thus, 

JDR 

we have that for v € H 2 (B R ), 

S(u,v)= \ fttV • ((ft • Vv)Vm - (Vm • Vv)ft) + £ 2 V • (huv)] dx - A AAi \ V • (VuV)dx 
J b r Jn 2 

where we have used the fact that h(x) = in the neighborhood of dB R . By the divergence theorem, the 
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right-hand side of the above equation becomes 

- J [((h ■ Vv)(Vw) + - ((Vw)+ • Vv)h) ■ v + k^h v uv] ds 

+A \ \((h ■ Vv)(Vw)_ - ((Vw)_ • Vv» ■ v + k\h v uv\ ds 
Js 

-^o J [((h ■ Vv)(Vm)+ - ((Vw)+ • Vv)/j) • v + fcf /i v wv] ds 

+A Ai J [((h ■ Vv)(Vw)_ - ((V«)_ • Vv)/z) • v + *f/i v wv"] - ^oA^i J ~f~ Vds 

-*JC 



- | - K{{Vu) + - A {Vu)_) ■ Vv + (*g - i £i)Mv 

— /<?«+ 



(A • Vv) I — - A x — I - h v ((Vu) + - /ii(V«)_) • Vv + (Jfc? - Ai*?)Av«v 

0V 0V ' 



-Vfifs 



(is 



-A0M1 f ^-T 
Js, dv 

= I [Av((V« + )r - Ao(Vb_) x ) • (Vv)/ - (*o - *)*i)Mv| 
-'So 

+/l f [^v((Va + )t - Ai(Vm_),) • (Vv), - (jfe j - A^hyuv] ds - AqAAi f -4~vds, (A.6) 

where we have used the transmission conditions ( 12.91 ) and ( 12.101 ) on 5;, i — 0, 1. Then we conclude that 
w' satisfies the equation d3.ll) in B# P| tlo and the equations ( 13.1D and ( 13.31 ) and that Lm' = du'_/dv\ dB . 
We extend m' outside of Br to be the solution of the Helmholtz equation AU + k^U = satisfying the 
Dirichlet boundary condition U = u'_ on 3Br and the radiation condition (I1.5I ). Denote the extension of 
m' by u' again. Then w' satisfies (13.11 ) in Qo> and we have 

r idu' du'_\_ r tdu' + du'_\_ 

S(u > v) = - J So 1^7 - Ao ^r s - Ao J Sl l^r - M -fo) vds - 

From this and dA.61 ) it follows that 



du' + du'_ 7 9i r \i 

~dv~ ~ ^"dv" = ° ~ A ° k ^ ho ^ vU \s + d/v So[( /i o)v((Vm + ), - Ao(Vu-),)\ on 5 , 

<9wl du'_ i i \ r ,i 

— i-^i— = (*f - **f)(fti)y«| Si +Div Sl [(A 1 ) v ((V M+ ) f -i 1 (V M _) J )J + 



^ du'_ 
dv 1 <9v 
The definition of u' gives 



A/li du + 



A[ dv 



on Si. 



u' - u'_ = -(ho) y 



1 du + 


du-\ 




\~dv 


~ ~dv) 


So 


I du + 


du-\ 




[dv 


~ ~dv) 


Si 



on So, 



on Si 



Since supp{h) c Br, we have ii(hMi) - u — w - U(h,AM) - u — u' on 3Br. This, together with the 
discussion just following (IA.5I) . implies that (u^ hAAl ^ - u s - m')q° - oCPollc'cSo) + H^iHcHSi) + l^il) as 
Pollc'(So) + H^illcHSi) + l^il ~* 0- We thus conclude that F is Frechet differentiable at (So, S\,X\) and 
the Frechet derivative F'\{s ,s^M)^^ h\,AA\) = u'^. The proof is then completed. □ 
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The following theorem can be shown similarly. 

Theorem A.3. Let Sq,S\ be C 2 and let u € Hj oc (^ 2 ) be the solution of the problem ft2.6\) - j[277\) , d2.9l ). 
rt2.77D -( |2".72D and ([73]) with T\ > 0. Then F is Frechet differentiate at (Sq,S \,ti) and the Frechet 
derivative F'|(S ,5i,ri) * s gi ven by F'\(SQ,Si,n){ho>h\, Ati) = u'^, where u'^ is the far-field pattern of 
u' e Hj oc (a ) H ff^Qi) H ^(Aj) w/i/c/i so/ves tfie problem (E2)-(E21> and (O wtfz 



/ <9w + 


du-\ 






~ ~dv) 


So 



fi = -0 

h = ( k o ~ ^i)(^o)v«| So + Div So [(h ) v ((Vu + ) t - Ao(Vu-) t )], 

* - Mf-fh 

U = (kj - ^!^)(/ii) vM | S] + Div Si [(/h) v ((Vk + ), - M<yuJ) t j\ - 
Remark A.4. By [9] again, u € H 2 oc {Q. Q ) f| # 2 (A) f| H 2 ^), so m' exists. 
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